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Abstract 

An example integrated flight/propulsion control system design is presented for the 
piloted longitudinal landing task with a modern, statically unstable fighter aircraft using 
the decentralized/hierarchical Design Method for Integrated Control Systems (DMICS) 
technique. The example design begins with an open— loop analysis on the scaled, linear 
integrated flight and propulsion system model. This leads to a partitioning of the 
integrated system model into separate engine and airframe subsystems, coupled only by 
their outputs. The airframe subsystem inputs are then converted to "generalized" controls. 
The airframe control is designed using a Linear Quadratic Regulator (LQR) based, explicit 
model following method for the "generalized" airframe. A control selector is designed to 
distribute the "generalized" control inputs over the actual airframe subsystem control 
inputs, which include the aerodynamic control surfaces and the engine— to— airframe 
interface variables that result from the partitioning step. The distribution of the airframe 
"generalized" control requirements over the interface variables (propulsive thrusts and 
moments) imposes performance requirements on the engine subsystem that are used to 
design a command— following LQ engine control. These imposed requirements from the 
high-level to the low-level system define the hierarchical nature of the design method. 
The performance and stability of the global plant with subcontrollers is evaluated and 
compared to that of the individual subsystems with their own independent controller. 


Controls Engineer. 
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Nomenclature 
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= thrust reverser and nozzle throat area (in ) 

= engine pressure ratio 
= expected value of [.] 

= total nozzle forces in the x and z direction (lbf) 

= open— loop plant transfer function matrix 
= controller dynamic gain matrix 
= engine fan and core speed (rpm) 

= mixing plane pressure (psi) 

= diagonal matrices for scaling plant input, output, and state vectors 

= high pressure turbine blade temperature (deg R) 

= total nozzle pitching moment (ft— lbs) 

= pilot input, selected airspeed (ft/sec) 


• 2 
V,V = aircraft forward acceleration (ft/sec ) and airspeed (ft/sec) 

WF36 = engine main fuel flow (lbm/hr) 

W ea ,r ea = airframe— to— engine interface matrices 

W M ,r M = engine— to— airframe interface matrices 
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= airframe and engine subsystem subscripts 
= commanded variable subscript 

= transpose of matrix or vector [•] 

= altitude (feet) 

= plant and scaled plant subscripts 

= body pitch acceleration (degs/sec ) 

= body pitch rate; state in (rads/sec); output in (degs/sec) 
= the Laplace variable 
= generalized variables superscript 
= body axis velocity (ft /sec) 

= input, state, and output vectors for the plant 


pTE = leading edge and trailing edge flap deflection (degs) 

t = pilot input, stick deflection (inches) 

«rv = thrust vectoring deflection (degs) 

A. = "i"th eigenvalue 

u = frequency in radians per second 

t = time constant (seconds) 

= pseudoinverse of matrix [•]. 


IFPC = Integrated Flight / Propulsion Control 
DMICS = Design Methods for Integrated Control Systems 
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Introduction 


Efforts to improve the maneuverability and enlarge the flight envelope of tactical 
aircraft via the use of propulsive moments and forces for flight control has led to an 
increase in coupling between the propulsion and airframe dynamics. This coupling has led 
to the inadequacy of the traditional approach of designing the flight and propulsion control 
systems separately. An integrated approach to flight /propulsion control system design is 
required to obtain an overall system which will yield improved performance as describe 
above, while reducing the pilot workload. 

In the early 1980’s, the U.S. Air Force initiated the Design Methods for Integrated 
Control Systems (DMICS) study with the objective of developing Integrated Flight/ 
Propulsion Control (IFPC) design methodologies for advanced tactical aircraft. The 
DMICS study resulted in two different approaches to the IFPC design problem: 

i) a decentralized, hierarchical approach which consists of partitioning the integrated 
system into subsystems and then designing separate controllers for each subsystem 
such that the high-level performance criterion are met [1]; 

ii) a centralized approach which consists of designing one global compensator at each 
operating point for the integrated system, using a Linear Quadratic Gaussian / Loop 
Transfer Recovery (LQG/LTR) design technique [2]. 

More recently, a joint NASA Lewis and NASA Ames IFPC research program [3] has been 
initiated with the objective of developing the technologies that will lead to a demonstrator 
aircraft with Supersonic, Short Take-Off and Vertical Landing (SSTOVL) capability in the 
1990’s. Among the goals of this program are the development of alternatives to the 
DMICS methodologies, improvements in overall system performance, and simplification of 
the control law synthesis and implementation. Prior to considering alternative techniques, 
the DMICS methodologies are being reviewed to develop an understanding of the control 
synthesis and evaluation procedures used and to identify the strengths and weaknesses of 
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these methodologies. Towards this goal, results are reported in this paper from an example 
application of the hierarchical, decentralized DMICS approach on a modern, statically 
unstable fighter aircraft for the piloted, longitudinal landing task. 

In the following, the decentralized, hierarchical DMICS procedure is first 
summarized and its application to the example problem presented in this paper is 
discussed. Then, the vehicle model used in the example study is described and the 
procedure for partitioning the integrated system is presented, along with a description of 
the resulting subsystems. The high-level airframe performance specifications and control 
design are then presented and the control performance is evaluated. Next, the generation 
of the low— level (engine) subsystem specifications from the airframe requirements are 
discussed. The engine performance specifications are presented along with the subsystem 
control design. A compensator to accommodate the influence of airframe outputs on the 
engine subsystem is also considered. Finally, the entire closed-loop system performance 
and stability characteristics are examined. The paper concludes with a discussion of the 
strengths and weaknesses of this decentralized, hierarchical approach to the design of 
integrated flight /propulsion control systems. 

Control Law Design Methodology 

The following is a discussion of the key features of the hierarchical, decentralized 
DMICS procedure, based on references [1] and [4]. Figure 1 describes the flow of 
information required for this design method and shows that the first requirements are a 
vehicle and mission definition, and a nonlinear simulation of the vehicle. Assuming that 
these are available, the procedure begins by obtaining an integrated airframe and engine 
linear model by perturbing the nonlinear simulation of the system*. This linear model is 
then scaled (normalized) to perform observability, controllability, and modal analyses. 


In reference [1] the integrated linear model was built from the individual subsystem linear models. 
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These analyses are used to establish modal groups of inputs, outputs, and elements of the 
state vector for partitioning the integrated system into subsystems. Scaling is important in 
making relative comparisons between variables, as will be shown later in the paper. The 
choice of proper scaling requires physical knowledge of both the airframe and propulsion 
systems and a detailed description of the nominal operating point about which the 
integrated linear model is obtained. 

After the variables are grouped into airframe and engine subsets, the system 
matrices are rearranged to order the variables into respective groups, resulting in the 
following matrix structure: 
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The subscripts "a" and "e" represent variable groups for the airframe and engine, 
respectively. More subsystems are possible, but only engine and airframe partitions are 
considered in this study 5 . It is desired to partition the system of equations in (1&2) such 
that the resulting subsystems are coupled only through their outputs (i.e., there is no direct 
coupling of the elements of the state vector). This desired structure is shown in Figure 2 
and corresponds to the following set of equations: 
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5 In reference [1], the propulsion subsystem was partitioned into inlet, gas generator and nozzle 
subsystems. 
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r and W are the interface matrices between subsystems and the nonsquare matrices V and 

a 

V e select the interface variables z & and z g from the available outputs y & and y . The input 
format in equation (3) shows that both the airframe physical controls (u ) and the 
engine— to— airframe interface variables (z^) are now considered as inputs to the airframe 
subsystem, while in equation (4), the airframe-to-engine interface variables (z ) are 
considered as disturbances to the propulsion subsystem. The above format is used to 
emphasize that the interaction from the airframe to the engine is not addressed directly in 
this methodology. 

The objective of partitioning is to select a minimum set of interface variables z and 

e 

\ such that the "errors" in representing the system in (1&2) using (3&4) are "small", in 
some sense. The selection of z & and z g are not unique. Any combination of variables that 
solves the model matching problem is satisfactory. One possible solution to this model 
matching problem is to minimize the square error in the matrix elements as determined by 
comparing the system described in (1&2) to the system description obtained by 
transforming the partitioned approximation of the system (3&4) back to the integrated 
structure (1&2). This approach requires that the dynamic characteristics of the system be 
relatively insensitive to small changes in the individual elements of the system matrices. 
This model matching problem requires the solution of 16 matrix equations in 12 unknown 
matrices, resulting in 4 matrix constraint equations, as described in the appendix. 

Once constrained to the structure shown in (3&4), the real problem in matching 
(1&2) is in the selection of a minimum number of variables and z ; that represent the 
interface between the subsystems. The interface represents the information flow between 
subsystems. If the effect of this information flow is not "small", there may be a 
requirement for a corresponding flow of information between the subsystem controllers in 



order to obtain good performance and stability properties. Therefore, the number of 
interface variables should be kept small to simplify the control design and implementation. 
Due to the way nonlinear aircraft simulations are currently developed, a set of outputs 
defining the physical airframe/propulsion interface can be obtained from the software 
interface between the engine and airframe nonlinear simulations. Once a set of interface 
variables is identified and a reasonable partition is obtained, the airframe subsystem (3) is 
taken as the high-level plant, while the engine is considered to be a low-level subsystem 
(4). This hierarchy is described by Figure 1, as the airframe subsystem becomes the 
mission level subsystem, and the engine becomes a function level subsystem. In Figure 1, 
the information flow between the Mission— Level and the Function— Level refers to the 
requirements and capabilities of the interface variable z ' in terms of its closed— loop 
z 

response . In the example presented here z is comprised of propulsive thrust and 
z 6 

e 

cmd 

moments. 

Now that the separate subsystems are defined by the partitioning step, the design of 

the high-level airframe controller can begin. First, the high-level airframe model is 

reformulated as a system with "generalized" control inputs [1]. There is one generalized 

control input for each degree of freedom for the rigid body aircraft. For the longitudinal 

model used in this study there are three degrees of freedom: forward and vertical 

displacements and a pitch rotation. These yield three generalized control inputs, £u, <5w, 

and 6q, for the forward, vertical, and pitch acceleration equations, respectively. The 

generalized control inputs are scaled such that a 100% deflection yields a specific value of 

acceleration, based on performance specifications, maximum control authority of the 

high-level system, and estimates of the low-level subsystem’s closed— loop performance. 

This defines the new input matrices B and D , which will yield a new "generalized" plant, 
* * * 

G (s)=D +C (si— A )B . Estimates of the nominal bandwidths of the generalized 

ci a a a' a ° 
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controls are made from the known dynamic characteristics of the actual physical controls of 
the high-level system and the estimated closed— loop response of the low-level systems. 
Referring back to Figure 1, the estimated capability of the closed — loop function level 
subsystems are used to define the available authority and bandwidth of the generalized 
controls. For example, u is mainly dependent on the available horizontal thrust, which is a 
capability that the engine provides. Thus, the bandwidth for the forward acceleration 
generalized actuator is defined by the estimated closed— loop horizontal thrust response. 

The regulator design can now begin on the generalized plant G a (s). The 
methodology does not dictate the technique to be used for the control design. In this paper 
the airframe control design is formulated as a stochastic tracking problem [5] and an 
explicit model-following Linear Quadratic Regulator [6] is designed for the airframe with 
generalized controls. The desired flight handling specifications are built into the model to 
be followed, called the Maneuver Command Generator (MCG) in reference [1], The basic 
LQR problem can be extended to provide for frequency shaping of the performance index, 
but this was not used in the example study for simplicity. Frequency shaping would affect 
the nominal bandwidth of the generalized actuators and thus it would also affect the 
performance requirements imposed on the engine subsystem. Once the LQR design yields 
an acceptable state feedback controller, it is approximated by output feedback, since not 
all the states are available as measurements. 

A control selector, T(s), is designed to distribute the generalized control inputs over 
the airframe inputs for the longitudinal case show below: 
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( 6 ) 


Figure 3 shows the resulting airframe control structure. This structure will be discussed in 
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detail later in the paper. Evaluation of the closed— loop system stability and performance 
completes the first iteration of the high-level control design. This controller design is now 
used to generate specifications for the engine subsystem. 

The engine subsystem specifications are derived by determining the level of 

uncertainty in the generalized inputs that can be tolerated for stability of the high-level 

system. This stability robustness bound is then translated from the generalized control 
— * 

inputs (u ) to both the high-level system inputs, which consist of the airframe physical 
controls (u a ) and the engine interface variables (z^). Reference [4] formulates this 
conversion as a model matching problem and presents the solution in the frequency 
domain. The net result is a set of nominal "actuator type" frequency response curves for 
the subsystem interface variables (z e ) and a set of bounds around the nominal response 
that the subsystem response must fall within to guarantee the stability of the integrated 
system with separate controllers. For good performance, the subsystem response must 
match the nominal closely. 

Now that the performance specifications for the engine subsystem are set, the engine 
compensator can be designed. Again, the methodology does not dictate a specific design 
approach. In this example, the subsystem design is similar to the high-level design in that 
a command— following LQ regulator is used. If the subcontroller design cannot meet the 
performance specifications, then the specifications need to be relaxed. This would require a 
bandwidth reduction for the generalized actuators, the redesign of the high-level regulator, 
the generation of a new set of subsystem specifications, and the redesign of the propulsion 
control. After the propulsion control design is complete, a disturbance rejection controller 
is considered to accommodate the influence of the airframe on the engine. This disturbance 
rejection controller was not required in this example problem because of the relatively weak 
airframe— to-engine coupling, as will be shown later. The resulting closed— loop engine 
subsystem control structure is shown in Figure 4. 
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After evaluating the closed— loop performance of the engine subsystem controller 
design, the entire integrated system model response with independent subcontrollers is 
evaluated. The performance of the integrated system is also compared to that of the 
separate individual subsystems. The resulting control structure for the integrated system 
is shown in Figure 5. 


Vehicle Model 


The vehicle model consists of an integrated airframe and propulsion state— space 
representation of a modern fighter aircraft powered by a turbofan engine and equipped with 
a 2D thrust vectoring nozzle. The vehicle dynamics are linearized at a flight condition 
representative of the Short Take — Off and Landing (STOL) approach to landing task 
(airspeed V^= 120 knots, flight path angle 7 0 = -3 deg.). The open— loop model is the same 
model as the one used in reference [7], with the addition of a few outputs. The vehicle 
model is defined as follows: 


x = A x +B u 
P P P p p 


= C x + D u 
P p p p 


where x p , u p , y p are the perturbed state, input, and output vectors as described below: 


x p = [u,w,q,0,h,N2,N25,P6,T41B] T 


(7) 


y p = [q,q^,7,7,V > V,h,F cx ,F ez ,T em> N2,EPR] T . 

The vector y p represents the outputs that were necessary to perform the design using this 
methodology. Not all the available outputs are shown and not all of the outputs shown are 
controlled variables. Additional outputs would be necessary in the actual implementation 
for rate and limit logic (rotor speeds, stall margin protection, etc.). The vehicle system 
matrices A p , B p , C p , and D p are listed in the Appendix. The open— loop vehicle 
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eigenvalues are 

X 1 =+0.073, A 2 3 = — 0.094±j0.23, A 4 =+1.07, A 5 =-1.47 (Airframe Modes) 

A g =—1.40, A 7 =— 3.57, A g =— 6.96, Ag=— 89.28 (Propulsion Modes) 

Note that the airframe is statically unstable with a highly unstable pitch mode. The 
open— loop plant transfer function matrix G^(s) was full, indicating that there is coupling 
between the airframe and the propulsion system. However, as discussed later, the 
airframe— to-engine coupling is weak for this particular model. 

Analysis of the control distribution matrix B p indicates that the leading and trailing 
edge flaps, 5 pLE and ^ FTE , are direct— lift devices that provide indirect control of the flight 
path angle. The fuel flow (WF36) thrust reverser area (A78) and nozzle throat area (A8) 
affect the engine dynamics, and also indirectly affect the airframe dynamics through the 
changes in the forces and moments. The thrust vectoring angle fl TV , while affecting the 
engine resultant thrust direction, is also the primary pitch control. 

System Partitioning 

Modal, controllability, observability, and steady— state gain analyses were performed 
to obtain the measures used to group the system variables. The V'th normalized column 
of the modal transformation matrix T. shows how the "i"th mode, C, contributes to the 
state vector x=T.<^. Similarly, a column of the control effectiveness matrix in modal 
coordinates, (T ^B)., shows how the corresponding input affects the modes, and a row of 
the output matrix in modal coordinates, (CT)., shows how the modes contribute to the 
corresponding output. The steady— state gain matrix is calculated after first stabilizing the 
system with an LQ control design using a heavy input cost. The unstable poles of the 
resulting closed— loop system were reflected about the imaginary axis and the stable poles 
remained fixed. The steady— state influence of the inputs on the outputs was also used to 
group the inputs and outputs. All of the measures described above are relative to the 
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nominal operating values about which the linear model was obtained. Xo allow the system 
variables to be compared on an equal basis, each input, output, and state element was 
scaled so that the "expected" linear range perturbation of the scaled variable would fall in 
the range of *1. This allows, for example, a change in flaps of 2 degrees to be directly 
compared to a 500 lbm/sec change in the fuel mass flow rate. The relative magnitudes and 
units of the variables are taken into account in the scaling. The resulting scaled system 
appears below: 


x = A x + B u 


P. P a 

S 3 


P « Pc 
8 8 


y = C x + D u 


p « p» 
8 8 


P P 
8 8 


( 8 ) 


where x p ^-S x x p , u p -S u u p , and y p =S y y p . The scaling matrices appear in the 


appendix. 


Tables 1 and 2 indicate the effect of scaling the modal measures used for grouping 
the state variables. The tables show the modal contribution to each state by displaying the 
normalized eigenvectors for each eigenvalue. For example, in Table 2 the fastest 
eigenvalue, A = -89.3 is most evident in the mixing plane pressure (P6) and A = +1.07 
contributes most strongly to q and 0 , the unstable pitch mode. In Table 1, the normalized 
"engine modes" have zero affect (to 3 decimal places) on the aircraft state elements 
u,w,q,0, and h. Comparing Table 1 to Table 2, the contribution of the "airframe modes" 
have little effect on the scaled engine state vector as compared to the effect on the the 
unsealed system. Consider the relative contribution of the complex airframe mode 
(— 0.094±i0.23) to fan speed (N2). This contribution drops from 0.582 to 0.06 after scaling 
the system. This shows that scaling is necessary for the direct comparison of these 
measures of modal contribution. Examining Table 2, it is clear that the ease of grouping 
the elements of the state vector results from the existence of decoupled modes. This 
indicates that there is a relatively small amount of coupling between the airframe and 
engine state elements in this model, which is also apparent from the scaled system matrix, 
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Controllability and observability analyses on the scaled system yield similar results. 


Table 1. Normalized Columns of the Eigenvectors of the Unsealed Open— Loop System 
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Table 2. Normalized Columns of the Eigenvectors of the Scaled Open-Loop System 
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During the open— loop analysis of the control input effectiveness, it became apparent 
that the leading and trailing edge flap deflections ^ FLE and ^ FTE belong with the airframe 
subsystem, and that fuel flow (WF36), thrust reverser port area (A78), and nozzle throat 
area (A8), belong with the engine subsystem. The placement of the thrust vectoring angle 
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(6p V )> was n °t so cosily determined. 6^ is the primary pitch control because it 
determines the direction of the thrust which generates the pitching moment. The physical 
engine-to-airframe interface is comprised of forces and moments. Thus, it was concluded 
that propulsive moments and forces are a reasonable set of variables to represent the 
physical engine-to-airframe interaction. Since S Ty determines the direction of the thrust, 
it was grouped as an engine control input. An alternative set of interface variables could 
been used (mass flows for example) but from a modelling viewpoint the optimal set of 
coupling variables for the partitioned linear models is the same set of variables used to 
represent the interface between the airframe and engine in the nonlinear simulations. 
From a simulation stand point, there is no other set of variables that will better represent 
the interface information. Thus, one possible set of interface variables for partitioning the 
linear system is predefined by the nonlinear simulation from which the linear system was 
obtained. 

Inlet distortion effects are one source of airframe— to-engine coupling. Although 7, 
q, and 0 contribute to this coupling, the primary effects are due airspeed (V) and altitude 
(h) for this example problem . Mach number and pressure can be used to schedule the 
engine operation, and they convey the same information as V and h. Therefore, the 
variations in V and h would be implicitly accounted for in the nonlinear portion of the 
propulsion control (the logic or schedule portion). For this particular model, the 
airframe-to-engine coupling is relatively small, but this will not be the case for high 
performance aircraft envisaged for the future. 

The partitioning study resulted in the following groups of variables: 
airframe state vector: x & = [u,w,q,0,h] engine state vector: x = [N2,N25,P6,T41B] T 
airframe inputs: \ = [^ fle ,^ fte ] T engine inputs: u = [WF,A78,A8,<* tv ] T 
airframe interface: z = [V, 7, q,0,h]T engine interface: z = [F ,F ,T 1^ 

a e L ex’ ez’ em J 

The key decision of this partitioning step was the placement of with the engine 
subsystem, as was previously described. Note that q, 7, V, N2, and EPR are not necessary 
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to describe the coupling between the engine and the airframe, z and z contain all the 
information necessary to describe the information flow between the subsystems, q, 7, and 
V are not needed, because the engine inlet conditions do not change as fast as these 
variables. N2 and EPR are not needed, because they are already represented in the thrust 
and moment variables. 

With the variables grouped, the system can be represented by (1&2) and then 
transformed to the matrix structure represented in (3&4). This structure defines the 
airframe partition and the high-level regulator design can begin on the system described by 
equation (3). 


Airframe Compensator Design and Evaluation 
The design specifications for the high-level airframe control system are as follows: 

(1) Track airspeed, flight path and pitch angles in a decoupled manner with zero 
steady— state error for step commands. The desired control band widths are 1 
rad/sec for the velocity loop and 5 rad/sec for both the pitch and flight path angle 
loops. 

(2) Avoid "excessive" control deflections and rates to prevent nonlinearities due to 
control deflection and/or rate limiting. 

(3) Maintain "adequate" stability margins in all control loops to guarantee stability in 
the presence of unmodelled dynamics and variation in model parameters. 

To achieve these goals, an LQR based, explicit model— following control system was 
designed. The model to be followed, called the Maneuver Command Generator (MCG) in 
reference [1], was designed to meet military specifications for Level I handling qualities 
requirements [8]. The inputs to the MCG, u mcg , are the pilot stick deflection (<5) s , and 

selected airspeed (V^); the MCG outputs, y mcg = [4>Q>0,7,7,V,V]^ cg , are the reference 
values for the airframe outputs to track. 
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The desired responses to pilot command inputs are listed in Table 3 in transfer 
function form. In Table 3, the desired velocity response represents a well-damped response 
with minimal overshoot and a rapid settling time with a 90% rise time to a step input 

q , 

4 rg = 5 secs. The selection of the pitch rate response /^ t is based on desired short 
period characteristics for Level I handling qualities. Also, open— loop analysis revealed that 
it would not be possible to control flight path angle independent of the pitch angle. 
Therefore, the flight path response is such as to lag the pitch attitude response with a time 
constant of 0.52 seconds. This corresponds to r^=0.52 in the "classical" airframe 
longitudinal dynamics [9]. It was desired to have t ^> 0.7, but after open-loop analysis of 
the flaps, it was discovered that this would not be possible. These dynamics shown in 
Table 3 can also be realized in a state space representation (A , B C D ) 

' meg’ meg’ meg’ meg' - 

Table 3. Desired AirFrame Response Transfer Functions 


cmd 


V 


s c 1 


Notation: { K(l/r)[C;w n ] s K(s+l/r)(s 2 +2Co> n s+u£) } 

0.04(3.13) q cmd 35.12(0.52) 1 = 0 1.92 

[ 0 . 89 ; 0 . 36] ' ^TT [0.89;2.24] ’ ~\t ~\t' (0.52) 


The next step in the design procedure of the high-level system is to restructure the 
plant inputs as generalized control inputs. For the longitudinal model there is one 
generalized control for each acceleration term: body axis forward acceleration (£u), body 
axis vertical acceleration (tfw), and pitch angular acceleration (tfq). For this example, the 

3|c 

generalized control input matrices B and D were determined based on maximum 

a a 

available accelerations in each direction using maximum deflection for the physical 
controls. This procedure was based on the discussion in reference [10]. Using u max = [15, 
15, 5000, 50, 100, 10] resulted in the following generalized control input matrices: 
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0 0 — 3.32e— 1 


1.49e-l 0 

0 1 
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0 3.76e— 2 
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0 

0 
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0 0 
0 0 

5.8e-3 

0 

* 

D = 

a 

— 8.77e— 3 1 . 05e— 2 
0 0 

0 
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0 0 

0 J 


— 1.46e— 1 -7.3e-3 

0 





0 

0 

0 



At this point in the design, generalized actuators would normally be included into 
the generalized airframe. Estimates of the bandwidths of the generalized actuators are 
made from the known bandwidths of the physical controls and the estimates of the 
closed— loop response of the engine subsystem. However, generalized actuators were not 
used in this example study to simplify the design. This implies an infinite bandwidth 
generalized actuator, which imposes an infinite bandwidth requirement on the response of 
the engine-to-airframe interface variables F /F cmd , F /F cmd , and T /T cmd . This is 
synonymous to designing a control for a plant without actuators and then verifying that 
the actuators do not affect the design. Although this overlooks an important step in the 
methodology, it simplified the design process while still allowing the methodology to be 
exercised, which is the primary purpose of this paper. 

The structure of the open— loop design plant needs to be defined. To generate the 
error vector for the regulator design, the airframe outputs are selected to be the same as 

the MCG outputs y & = [q,q,0,7,7,V,V] . Lead information is obtained from the MCG 

outputs q, 7, and V, but these will not be included as airframe outputs in the output 
feedback control design. Next, the specifications call for zero steady— state error for step 
commands, requiring integral action on 3 loops: pitch attitude (0), flight path angle (7), 
and airspeed (V). Also, it was desired to design the controller for a pilot stick bandwidth 
of 10 rad/sec and a pilot airspeed select bandwidth of 1 rad/sec. The pilot input 
bandwidths were modeled using first-order lags of white noise as shown below: 

Kr 10t >6 “ d v «r -i v . d + i-V ( 10 ) 

where rj ^ and u v are zero mean, Gaussian, white noise processes. Combining the 
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generalized vehicle dynamics, the model dynamics, the integral errors and the pilot input 
models results in the following open— loop design structure: 
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( 12 ) 


The superscript "s" indicates a submatrix, since the integral errors on [0,7,V] are a subset 
of the error vector e = y mcg - y a = [e^, e q , e^, e^, e^, e* , e y ] T . The zeros are zero matrices 
of appropriate size. The subscript "a" describes the generalized airframe. A,, and B 

filter filter 

describe the state space model of the pilot inputs u mc g described in equation (10). 

Minimizing the following performance index results in an optimal state feedback 
compensator K^: 

T 

J = E{ 1 im - / ( e T Q e + u T Ru) dt } (13) 

T-»® T o ^ 

for a system of the form 


x = Ax+Bu+r£ y = C x + Du (14) 

ooooo 00000 ' ' 

where £ is a zero mean, Gaussian white noise process with identity intensity, r is the noise 
distribution matrix, y Q are the output measurements, and Q (positive semidefinite) and R 
(positive definite) are symmetric weighting matrices. The subscript "o" indicates the 
open— loop design plant described in (11&12). The optimal state feedback solution yields a 
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set of feedback gains corresponding to the design state vector as shown below: 

i x3 3x2 

: r ,K ] . 

Je’ u J 


u =K -x = [K ,K , 

» ▼ « L a 7 meg 
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The gains on the pilot inputs are feedforward, while the remaining state feedback gains 

p 

K a , K mcgJ and K ^ are transformed to output feedback gains resulting in a set of gains on 
the system outputs as shown below: 


u = K • 

a y 


r y 

* meg 
V 

3x7 3x4 3x3 

y 1 

_mcg 

v 

^ e rr 

= [K ,K ,K ]. 
y y yrj 

meg err J e 

y err 

p/«j 

l*M 


(16) 


The gain matrix K is dimension 3x4 because it does not include the outputs q, 7 , and V 

y err 

as stated previously. The state feedback conversion to output feedback was accomplished 
using matrix algebra and the resulting response matched the state feedback response. The 
resulting control structure is shown in Figure 3, which shows that this control law design 
results in a multivariable, proportional plus integral (PI) controller with additional 


feedforward gain matrices K and K . The integral gains K are only on the 7 , 0 , 

U P y mcg y Je 

and V loops, while K is the matrix of proportional gains. In evaluating the performance 

y err 

of the control design, the generalized control inputs were checked for rate and magnitude 
limits for step response inputs to 5 gt = 1 inch and V gel = 20 ft/s. The closed— loop 
bandwidths were acceptable for all loops. 

To this point, the controller has been designed for the generalized plant with 
generalized control inputs. After acceptable output feedback control performance has been 
obtained on the generalized plant, the next step in the methodology is to design a control 
selector which transforms the required generalized controls to the available physical 
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controls, as in equation (6). The control selector can be complex (including dynamics, 
reconfiguration and fault accommodation logic, etc.), but to simplify the design, a constant 
gain control selector was used in this study. The idea is to transform the generalized 

t 

control inputs to the physical airframe inputs and the subsystem interface variables as 
shown below: 


* * u 

Bu = B - 
a a a z 


A simple solution is to minimize the square error using a pseudoinverse, which results in 
the following: 


U a Jt *_* 

- = B *-B u = Tn 


Scaling the variables before using the pseudoinverse yields an even distribution of the 
relative residual errors and a control selector design that distributes the generalized control 
requirements uniformly over the available control effectors. During the design, it was 
discovered that this particular aircraft model had been linearized at an operating point 
where the propulsive lift port was closed. Thus, propulsive lift is limited to the thrust 
vectoring nozzle and the total propulsive lift will be less than what should be available 
from this aircraft. For this reason, the control selector was then modified such that it does 
not request any lift from the propulsion system (F^™^=0) by setting the corresponding 
elements of T to zero. This places an additional lift burden on the airframe, which results 
in large flap deflections for certain maneuvers. 


The closed loop stability and performance of the high-level controller with the 
control selector were evaluated on the nongeneralized plant and found to closely match 
that of the generalized plant. The resulting frequency responses of V/V and q/6 are 
compared to the desired MCG frequency responses in Figure 6 and 7. Figure 6 shows that 
the frequency response of the aircraft airspeed to the pilot velocity select command closely 


20 



follows the desired MCG response out to the bandwidth of interest. Figure 7 shows that 
the pitch rate response closely matches the desired MCG frequency response, until about 2 
rad/sec where the magnitude response runs about 2—3 db larger than the desired MCG 
response. This does not degrade the time response, as will be shown later. 

Subsystem Specification Generation 

The generalized actuators describe a set of nominal bandwidth requirements for the 
generalized controls. This bandwidth requirement can be translated from the generalized 
controls to the physical airframe controls (u^) and the interface variables (ij through the 
control selector (T). The level of uncertainty that the high-level system can tolerate for 
stability now needs to be calculated. Lehtomaki et al. [11] have shown that for 
multi— input multi-output systems, the minimum singular value of the return difference 
matrix a[I+KG] is a reliable measure of closed— loop system stability robustness to 
unstructured uncertainties occurring at the plant input. _ This uncertainty on the 
generalized inputs can be translated to the physical airframe and the airframe— to— engine 
interface variables as described in [4]. The result is a set of nominal bandwidth 
requirements for F /F^™ d , F /F cmd , and T /T cmd , and a tolerance around each 
nominal response denoting the response requirement for stability of the high-level system. 

This example study was simplified by not using generalized actuators, which implies 
an infinite bandwidth, which was discussed previously. Rather than calculating a stability 
bound on the subsystem response as in [4], the change in robustness due to the addition of 
the lower level system into the integrated system will be shown. 

Propulsion Compensator Design and Evaluation 
The gross thrust level of a turbofan engine can be controlled by using feedback control on 
the fan or compressor rotor speed and the engine pressure ratio [12], Based on engine 
thrust performance requirements and open— loop engine control effectiveness, the desired 
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control bandwidths were chosen as 5 rad/sec for the N2 loop and 10 rads/sec for the EPR 
loop. The control structure used in this study assumes the propulsion system receives 
setpoint information for N2 and EPR from an engine schedule as shown in Figure 4. The 
interface variables F and T (F cmd =0 as discussed earlier) are included in the 
closed— loop design, but integral control of the interface variables is the responsibility of the 
high-level system. For example, if there is a high-level velocity error, the error integral 
will accumulate and increase the request for propulsive thrust from the engine subsystem. 
The high-level system drives the request for thrust and there is no need for integral control 
of the net moments or forces in the engine subsystem. A feedforward gain matrix was used 
to compensate for steady— state errors at the nominal operating point. The specifications 
for F gx and T gm were taken from the thrust vectoring actuator bandwidth. 

Table 4. Desired Engine Response Transfer Functions 

Notation: { K(l/r)[C;w n ] = K(s+l/r)(s 2 +2<a> n s+w 2 ) } 

N2 _ 5.0 , EPR 10.0 . 15.0 . ^ em _ 15.0 

N2 c m d< 5 -°)' (10-°)’ C d ( 15 °)’ 

The F^ and responses were only limited by this actuator bandwidth due to the direct 
feed-through terms in the engine subsystem "D" matrix. 

The engine control design was formulated as a command following problem. This is 
similar to the explicit model-following approach used in the high-level design, except there 
is no model to be followed, so the commands are directly tracked. Since zero steady— state 
error for step inputs was desired for the N2 and EPR loops, two integral errors were 
appended to the system equations. Actuator dynamics for fuel flow (WF36), thrust 
reverser area (A78), nozzle throat area (A8), and thrust vectoring (£py), were incorporated 
into the design plant model. The engine actuator dynamics were approximated with four 
first-order filters with time constants of 0.1, 0.033, 0.033, and 0.033 seconds [12], which are 
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realized by the state space matrices A and B . Models of the commanded input 

act act 

signals, u cmd = [N2 cmd , EPR cmd , F^ d , T^ d ]^ were approximated by four first-order 
filters with time constants of 0.1, 0.067, 0.067, and 0.067 seconds, respectively, which are 
realized in state space description by A fil and B fi j. The resulting open— loop design 
structure is as follows: 

e 
s 
e 

act 





where u e = [WF36, A78, A8, <5 Ty ] T and y err = errors on [N2, EPR, F ex , T em ] T . Again, the 
superscript "s" signifies a submatrix since only the integral errors of N2 and EPR are used, 
and these are a subset of the output errors = y cmd ~ 7 eng - Minimization of a quadratic 
performance measure for the system given by (19&20) by using the well known LQR 
solution yields a full state feedback controller. This controller was used to complete this 
study because all of the engine states can be measured except for T41B, which may be 
estimated using P6 and T7, the measurable nozzle inlet temperature. The full state 
feedback engine control was used to complete this study. The resulting control structure is 
shown in Figure 4 and equation (21) as 
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The frequency response plots for N2/N2 cmd , EPR/EPR cmd , F /F cmd , and 

cmd 

^em^em 3X6 s ^ own Figures 8 and 9. Note that all responses meet the specifications 
stated in Table 4. The EPR bandwidth 6 is actually larger than desired and an additional 
filter may be required to increase the rolloff at higher frequencies. This wide bandwidth for 
the EPR response was required in order to obtain the desired response for the fan speed due 
to the coupling between EPR and N2. The propulsive thrusts and moment responses have 
large contributions from the direct feed-through terms ("D" matrix) that are related 
directly to the thrust vectoring angle. Therefore, the responses of the propulsive forces and 
moment are only limited by the thrust vectoring actuators. This is why the responses for 
F gx and T gm roll off after a frequency of 30 rad/sec. 

With the engine regulator designed, attention can be focused on the accommodation 
of the airframe— to-engine coupling. By cross feeding the airframe— to— engine coupling 
variables (V,7,q, 0 ; and h), a "measurable" disturbance accommodation controller can be 
constructed as shown in Figure 4. After investigating the effects of the airframe— to— engine 
coupling, it was determined that this additional control was not necessary for this 
longitudinal model because this coupling was weak. To show this weak coupling the 
magnitude of the frequency responses for forward thrust (F ) to the engine inputs 
(WF36,A8,A78,V,7,0,h) are shown in Figure 10. Note that the largest effects due to h and 
V are small, relative to the effect of A78 and A8. The altitude (h) has the largest effect, 
but it would normally be part of the nonlinear engine control logic (i.e., engine operating 
point as a function of altitude). Also, this operating point is h = 100 feet, so the expected 
linear range for h is less than 100 feet. Thus, the disturbance accommodation filter was not 
required for this example. 


Here, we use the accepted measure for bandwidth as the —3 dB amplitude point or —45 degrees phase, 
whichever occurs at the lower frequency. 
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Global Compensated System Evaluation 

The separate subsystem controllers were used to control the integrated, global plant 
as shown in Figure 5. The closed— loop eigenvalues of the integrated system compare well 
with the eigenvalues of the individual subsystems. A comparison of the frequency 
responses of the pitch rate and airspeed to pilot stick and pilot velocity select are shown in 
Figures 6 and 7. Figure 11 shows the time response of the airspeed, pitch rate and pitch 
angle for a doublet pilot stick input. Note that the integrated system closely matches the 
desired response of the MCG. Figure 12 compares the time response of the integrated 
system to that of the airframe alone and to the MCG, for a step input to the pilot velocity 
select command. Note that the response to V se] is reasonably decoupled from the flight 
path angle. Also, the state of the engine is not largely affected by this step command. The 
net result is a stable system that tracks the Maneuver Command Generator responses to 
the pilot inputs. And finally, Figure 13 compares the singular values of the return 
difference matrix I+K(s)G(s), plotted versus frequency for the airframe and the integrated 
system. In this figure there are 6 plots, but only 3 appear because the others match so 
closely. This figure shows no increase in uncertainty, due to the addition of the propulsion 
subsystem dynamics (i.e. the additional thrust dynamics does not drastically affect the 
airframe performance). 


Conclusion 

The objective of this study was to gain an understanding of the Decentralized, 
Hierarchical Integrated Flight/Propulsion Control design methodology that was developed 
under the Air Force— sponsored Design Methods for Integrated Control Systems program. 
This goal was achieved by exercising the design methodology for the linear point design 
presented in this paper. The control system design was performed for the piloted 
longitudinal landing task of a modern, statically unstable fighter aircraft, powered by a two 
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spool turbofan engine and equipped with a 2D thrust vectoring nozzle. The results of this 

study lead to the following summary of the strengths and weaknesses of the linear design 

portion of this methodology: 

Strengths 

1) This hierarchical approach results in lower order subsystems with easy to implement 
subcontrollers, compared to a single, global compensator. This approach allows 
system specialists to work in their area of expertise. Also, the selection of design 
criteria for the separate partitions is simplified, when compared to the selection of 
design criteria for a global LQR design approach. 

2) The use of generalized controls is unique and relieves the pilot from the task of 
distributing moment and thrust requirements among the available effectors. The 
use of generalized controls may also reduce the number of gains required to be 
scheduled when the control is in full flight envelope operation. 

Weaknesses 

1) The procedure does not consider the airframe— to— engine coupling in the linear 
portion of the design. Therefore, the methodology may require modifications to 
ensure the stability of the integrated systems on future V/STOL applications with 
stronger airframe-to-engine coupling. 

2) Estimates of the propulsive moments and forces were assumed to be available for 
the engine feedback control system. The quality of these estimates is untested and 
needs to be evaluated in terms of accuracy and stability when implemented. 

3) The method requires many steps and could lead to a highly iterative design if the 
specifications for the low-level subsystem are not satisfied. 

4) The design of the control selector was simplified to a constant gain matrix in the 
example discussed in this paper. The design of a general control selector capable of 
handing the entire flight envelope would be a more difficult task. 
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A discussion of the current IFPC research program that will evaluate this methodology for 
a more strongly coupled system over a wider flight envelope appears in reference [3]. This 
program will also evaluate the use of thrust estimation for propulsion control. 


Appendix 

Derivation of Linear System Matrix Partitioning 

The following describes the derivation of the solution to the partitioning problem 
using the minimization of the least square matrix element error criterion. The 
transformation from (3&4) to (1&2) results in 16 equations in 16 unknowns and can be 
solved with simple matrix algebra. The mapping is one— to— one (assuming the inverses 
below exist), and offers the unique solution shown below: 


M = (I-W V W V ) -1 and 

N ea a ae e' 


N = (I-W V W V ) 

' ae e ea a' 


-1 


Note that: NW V = W V M and MW V = W V N 

ae e ae e ea a ea a 
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NW V D = W VMD = W VD„ 
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MW VD = W V ND = W VD 
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A 11= A+r V MW VC = a +r vc„ 

11 a aee eaaa a ac e 21 


a = r VMC 

12 ac c e 


a = r v NC 

21 ca a a 


A = A +r V NW VC = a +r vc„ 

22 e ea a ae e e e ea a 12 


B.^B+r V MW VD =B+r V D 01 

11 a ae e ea a a a ae e 21 


B = r VMD 

12 ac e c 


B = r V ND 

21 ca a a 


B„= B +r V NW V D = B +r VD„ 

12 e ca a ac e c c ca a 12 


(A.l) 

The inverse of this mapping results in 16 equations in 12 unknowns. Eight of the equations 
can be combined to yield 4 constraint equations. The defining equations for the mapping 
from (1&2) to (3&4) are shown below along with the least square solution to the constraint 
equations: 


A = A 

a 11 


r vc„ 

ae e 21 


c = c -w vc 

a 11 ae e 21 

a = A 09 -r v c 10 

e 2.2 ea a 12 
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B = B — r V D 

a 11 ae e 21 

D = D— W V D 

a 11 u e 21 

B = B-.-T V D 

e 22 ea a 12 

D e = D„„-W VD 


A 12 = r ae V e C 22 


22 "ea a 12 

Least Square Solution of 4 Constraints 
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(A.2) 


B 12= r ae V eD 22 ^ = [B^ AJ.{V [D^ CJ} 
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a = r v c 1t 

21 ea a 11 


b 91 = r v d 

21 ea a 11 


c = W V C.. D 01 = W V D,. 

21 ea a 11 21 ea a 11 


w « = I D 12 C 1 2 MV [D 22 C 22 ]> 

F ca = f B 21 A 2 J ’ l V af D ll C ll^ 

W « = [ D 2i C 2i]-{ v .[ D n c n]} # < A - 3 ) 


# 


The solutions in (A.3) can be used to solve (A.2) for A ,B ,C ,D . A ,B ,C ,and D . These 
solutions are dependent on the C„ and D„ matrices, and on the set of outputs selected to 
represent the coupling interface, defined by V and V . Thus, once the interface variables 
are selected, the structure in (1&2) is transformed into the structure in (3&4) using (A.2, 
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A.3) above. The accuracy of this transformation is checked by transforming the resulting 
structure in (3&4) back to the structure in (1&2) using (A.l). By comparing the results of 
this retransformation to the original matrices the individual matrix element errors can be 
identified. If the errors are "reasonably small" and the system eigenvalues and 
eigenvectors are insensitive to these errors, then the transformation is successful and the 
structure in (3&4) is a valid representation of the dynamics system in (1&2). This scheme 
worked reasonably well for the example presented here, but, in general, a more detailed 
analysis will be required. Current research in optimization is applicable to this model 
matching problem and should result in a more general technique. 


Integrated System Matrices Required for Partitioning 

6 “P” 18 “ = « fie ,« fte ,WF36,A78,AM tv 

8 outputs y = V,7,q,0,h,F F T 

ex ez em 

9th order state vector x = u,w,q,0,h N2,N25,P6,T41B 


A p = COLUMNS 1 THRU 6 

— 5.8927e— 02 1.0675e-Ol -3.8598e+01 -3.1839e+01 1.4098e-02 3.1445e-04 
— 2.6588e-01 -2.6652e-01 1.9481e+02 -4.5989e+00 5.1963e-04 -1.5785e-05 
— 1.5414e-03 7.8060e-03 -1.9486e-01 -^1.8182e-04 2.5644e-05 9.4632e-07 
0.0000e-01 0.0000e-01 1.0000e+00 0.0000e-01 0.0000e-01 0.0000e-01 
1.4275e-01 -9.8976e-01 0.0000e-01 2.0060e+02 0.0000e-01 0.0000e-01 
7.7818e-01 1.5424e-01 0.0000e-01 0.0000e-01 -8.4851e-02 -4.1914e+00 
1.5179e-01 3.0085e-02 0.0000e-01 0.0000e-01 -1.6551e-02 4.2632e-01 
7.9337e-01 1.5725e-01 0.0000e-01 0.0000e-01 -3.5024e-01 2.2953e-01 
— 1.0053e— 01 — 1.9925e— 02 0.0000e-01 0.0000e-01 1.0962e-02 3.7401e-02 

COLUMNS 7 THRU 9 
2.5989e-04 3.8186e-02 2.2508e-03 
— 2.1058e— 06 1.8261e-04 -2.9567e-06 
3.7441e— 07 3.6678e-05 2.6732e-06 
0.0000e-01 0.0000e-01 0.0000e-01 
0.0000e-01 0.0000e-01 0.0000e-01 
6.0216e+00 — 3.4337e+02 1.1603e+01 
— 5.7070e+00 2.7163e+01 1.0396e+01 
1.1548e-01 -9.0238e+01 8.4760e-01 
— 1.0362e-01 — 7.9538e+00 -1.0682e+00 
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Bp= 3.4680e— 02 — 4.9602e-02 3.4360e-05 -2.0546e-01 6.9121e-02 -4.1834e-04 

6.9252e— 02 — 1.4550e-01 1.2340e-08 -2.9360e-04 7.1041e-O5 -5.4520e-01 
— 8.1002e— 03 7.1325e-04 5.5072e-08 1.0677e-04 -8.3833e-05 -7.9733e-02 
O.OOOOe— 01 O.OOOOe— 01 0.0000e-01 0.0000e-01 0.0000e-01 0.0000e-01 
O.OOOOe— 01 0.0000e-01 0.0000e-01 0.0000e-01 0.0000e-01 0.0000e-01 
O.OOOOe— 01 0.0000e-01 1.4694e-01 0.0000e-01 0.0000e-01 0.0000e-01 
O.OOOOe— 01 0.0000e-01 5.3656e-02 0.0000e-01 0.0000e-01 O.OOOOe-Ol 
0.0000e-01 O.OOOOe-Ol 1.8127e-02 -4.3021e+01 -2.5835e+01 O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe-Ol 1.6430e-01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 

C p = COLUMNS 1 THRU 6 

9.7971e-01 1.9412e-01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 
5.8883e— 02 -2.7958e-01 O.OOOOe-Ol 5.7296e+01 O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe-Ol O.OOOOe— 01 5.7296e+01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol 5.7296e+01 O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe-Ol O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol 1.0000e+00 O.OOOOe-Ol 
1.2147e+00 2.4076e-01 O.OOOOe-Ol O.OOOOe-Ol 1.2048e+01 3.4557e-01 
— 2.7442e— 03 ^5.4390e-04 O.OOOOe-Ol O.OOOOe-Ol -4.1640e-02 -7.4322e-04 
2.7487e— 01 5.4480e-02 O.OOOOe-Ol O.OOOOe-Ol 2.2269e+00 8.1022e-02 

COLUMNS 7 THRU 9 
O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol 
O.OOOOe-Ol O.OOOOe— 01 O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol 
2.3118e-01 3.1521e+01 1.9262e+00 
— 8.3705e— 04 -1.3301e-01 -7.5608e-03 
2.8577e— 02 2.4726e+00 1.9388e-01 

Dp= O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 

O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol O.OOOOe-Ol 
O.OOOOe-Ol O.OOOOe— 01 2.9552e-02 -1.7079e+02 5.6842e+01 O.OOOOe-Ol 
O.OOOOe— 01 O.OOOOe— 01 -9.3010e-05 -2.8600e-01 1.4822e-01 -4.5932e+02 
O.OOOOe-Ol O.OOOOe-Ol 4.7079e-03 9.6000e+00 -7.8168e+00 -9.3011e+03 


Scaling Matrices used for partitioning 

S u =diag([10 10 5000 50 100 20]), W W)’ WF36 ’ A78 > A8 - ^TV 

S =diag([20, 4.5, 4.5, 4.0, 20, 4000,400,10000]), V, 7 ,q,<?,h,F ,F ,T 

j ^ ex ez em 

S x = diag([20, 25, 0.0785, 0.070, 100, 573, 727, 30, 50]), u,w,q,0,h,N2,N25,P6,T41B 
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Airframe Partition System Matrices 


5 inputs 
5 outputs 

5th order state vector 


U a ^FLE’^FTE ,F ex ,F ez ,T em 


z a = V,7,q,0,h 
\= u,w,q,0,h 


A a = -6.0381e-02 1.0646e-01 -3.8598e+01 -3.1839e+01 -3.2590e-04 

—2.6589e— 01 — 2.6652e— 01 1.9481e+02 -4.5989e+00 4.2667e-04 
— 1.5445e-03 7.8053e-O3 -1.9486e-01 -4.8182e-04 -8.7528e-07 
O.OOOOe+OO O.OOOOe+OO 1.0000e+00 0.0000e+00 O.OOOOe+OO 
1.4275e— 01 -9.8976e-01 0.0000e+00 2.0060e+02 -9.9000e-05 


[B T ] = 3.4680e-02 -4.9602e-02 1.2026e-03 4.2499e-04 -2.0943e-05 

L a ae J 

6.9252e-02 — 1.4550e-01 3.0735e-06 4.9203e-04 3.4319e-05 

— 8.1002e— 03 7.1325e-04 6.1723e-08 ^1.4099e-05 1.0750e-05 

O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 0.0000e+00 

O.OOOOe+OO O.OOOOe+OO 0.0000e+00 0.0000e+00 O.OOOOe+OO 


C = — 8.8495e-02 4.4721e-01 -1.1165e+01 -2.7606e-02 -5.0150e-05 

a 

O.OOOOe+OO 0.0000e+00 5.7296e+01 0.0000e+00 O.OOOOe+OO 
O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 5.7296e+01 O.OOOOe+OO 
7.0783e— 02 8.0783e-02 5.5794^-01 -5.8901e-01 -1.3848e-04 
5.8883e— 02 — 2.7958e— 01 O.OOOOe+OO 5.7296e+01 O.OOOOe+OO 
— 1.1077e-01 5.2564&-02 1.6706e-03 -3.2086e+01 -2.3647e-04 
9.7971e-01 1.94l2e-01 O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 


[D W ]= -4.6411e— 01 4.0866e-02 

a ae J 

0.0000e+00 O.OOOOe+OO 
O.OOOOe+OO O.OOOOe+OO 
— 1.7320e— 02 3.7759e-02 
O.OOOOe+OO O.OOOOe+OO 
4.7420e— 02 -7.6840e-02 
O.OOOOe+OO O.OOOOe+OO 


3.5365e— 06 -2.5267e-03 6.1594e-04 

O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 
O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 
6.9951e-05 -1.1254e-04 -1.0828&-05 
O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 
1.1788e-03 5.1188e-04 -1.3856e-05 
O.OOOOe+OO O.OOOOe+OO O.OOOOe+OO 
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Propulsion Paxtition System Variables 


4 inputs u = WF36,A78,A8,&_ 

c 1 V 

5 coupling inputs z = V,7,q,0,h 

3 outputs z = F ,F ,T 

e ex ez cm 

4th order state vector x = N2.N25.P6.T41b 


A e = — 4.1914e+00 6.0216e+00 -3.4337e+02 1.1603e+01 

4.2632e— 01 -5.7070e+00 2.7163e+01 1.0396e+01 
2.2953e— 01 1.1548e-01 -9.0238e+01 8.4760e-01 
3.7401e— 02 -1.0362e-01 -7.9538e+00 -1.0682e+00 

B c = 1.4694e— 01 0.0000e-01 0.0000e-01 0.0000e-01 

5.3656e— 02 0.0000e-01 0.0000e-01 0.0000e-01 
1.8127e— 02 — 4.3021e+01 -2.5835e+01 0.0000e-01 
1.6430e— 01 O.OOOOe— 01 0.0000e-01 0.0000e-01 

C e = 3.4557e— 01 2.3118e-01 3.1521e+01 1.9262e+00 

— 7,4322e— 04 -8.3705e-04 -1.3301e-01 -7.5608e-03 
8.1022e— 02 2.8577e-02 2.4726e+00 1.9388e-01 

D e = 2.9552e— 02 -1.7079e+02 5.6842e+01 0.0000e-01 

— 9.301 Oe-05 -2.8600e-01 1.4822e-01 -4.5932e+02 - 
4.7079e— 03 9.6000e+00 -7.8168e+00 -9.3011e+03 

r ca = 7-9431e-01 — 1.7582e— 04 0.0000e-01 1.7582e-04 -8.4851e-02 

1.5494e— 01 -3.0177e-05 0.0000e-01 3.0177e-05 -1.6551e-02 
8.0981e— 01 — 1.7668e-04 0.0000e-01 1.7668e-04 -3.5024e-01 
— 1.0261e— 01 2.0460e-05 0.0000e-01 -2.0460e-05 1.0962e-02 

W ca = 1.2399e+00 -2.7131e-04 0.0000e-01 2.7131e-04 1.2048e+01 

— 2.8011e— 03 5.6481e-07 0.0000&-01 -5.6481e-07 -4.1640e-02 
2.8057e-01 — 5.7632e-05 0.0000e-01 5.7632e-05 2.2269e+00 


Acknowledgement 

The authors would like to thank Mr. Kevin Madden and Dr. P.D. Shaw, Flight 
Control Research Dept., Northrop Corp., Hawthorne, CA., for having provided the 
linearized aircraft model used in this study. 


32 



References 


[1] Shaw,P.D., Rock, S.M., and Fisk, W.S., "Design Methods for Integrated Control 
Systems", AFWAL-TR-88-2061, Wright Patterson AFB, Dayton OH, June 1988. 

[2] Smith, K.L., "Design Methods for Integrated Control Systems", 

AFW AL— T R— 8 60— 2 103, Wright Patterson AFB, Dayton, OH, Dec. 1986. 

[3] Mihaloew, J., "Flight Propulsion Control Integration for V/STOL Aircraft", NASA 
TM— 100226, Dec. 1987. 

[4] Rock, S.M., Emami— Naeini, A., and Anex, R.P., "Propulsion Control Specifications 
in Integrated Flight/Propulsion Control Systems", AIAA-88-3236, 24th Joint 
Propulsion Conference, July 11, 1988. 

[5] Kwakernaak, H., Sivan, R., " Linear Optimal Control Systems ". Wiley Interscience, 
1972, p. 253-265. 

[6] "Multi-Input Multi-Output System Control for Experimental Aircraft", 
NASA— CR— 181197, Purdue University School of Aeronautics and Astronautics 
Report N37-26915. 

[7] Garg, S., Mattern, D., and Bullard, R.E., "Integrated Flight/Propulsion Control 
System Design Based on a Centralized Approach", AIAA Guidance, Navigation and 
Control Conference, AIAA 89-3520, Boston MA, August 1989. 

[8] "Military Specification - Flying Qualities of Piloted Airplanes", MIL-F-8785C, 
USAF, Wright Patterson AFB, OH, Nov., 1980. 

[9] McRuer,D., Ashkenas,I., and Graham, D., " Aircraft Dynamics and Automatic 
Control", Princeton University Press, Princeton, NJ, 1973. 

[10] US/Canada STOVL Technology Program, Monthly Technical Narrative No. 15, 
December 1988, by General Electric for NASA Contract NAS3-25193. 

[11] Lehtomaki, N.A., Sandell Jr., N.R., and Athans, M., "Robustness Results in 
LQG— Based Multivariable Control Design", IEEE Transactions on Automatic 
Control, Vol. AC-26, Feb. 1981, pp. 75-92. 

[12] Systems Control, Inc., "F100 Multivariable Control Synthesis Program", 

AFAPL— TR— 77— 35, June 1977. 6 


33 



VEHICLE 0/EHICLE 

DEFINITION CAPABILITIES 



LD 

< 


OJ 

o 


o 

o 


(JJ 

00 

a 


o 

o 


LU 

OS 


Z LU 

o as 

oo =3 
oo o 




34 


FIGURE 2. - BLOCK DIAGRAM FOR SYSTEM PARTITIONING STRUCTURE. 











FIGURE 3. - HIGH LEVEL AIRFRAME CONTROL STRUCTURE. 
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FIGURE - LOW LEVEL ENGINE CONTROL STRUCTURE. 



FIGURE 5. - DECENTRALIZED, HIERARCHICAL IFPC CONTROL STRUCTURE. 
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FIGURE 10. - A I RFRAME-TO-ENGI NE COUPLING. 
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FIGURE 11. - AIRFRAME RESPONSE TO PILOT STICK DOUBLET. 
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FIGURE 12. - INTEGRATED SYSTEM RESPONSE TO V SEL STEP = 20 FT/S, 
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propulsive system 
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FIGURE 13. - CHANGE IN SINGULAR VALUES DUE TO THE 
ADDITION OF THE ENGINE SUBSYSTEM. 
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